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We describe the contributions of virtual corrections and soft gluon emission to the inclusive Higgs 
boson production cross section pp — > H + X computed at next-to-next-to-leading order in the heavy 
top quark limit. We also discuss estimates of the leading non-soft corrections. 



I. INTRODUCTION 

The Standard Model is almost thirty-five years old, and its essential goal, to describe the electro- weak interac- 
tions as a spontaneously broken SU(2)l ®U(1)y gauge symmetry has been spectacularly confirmed. However, 
the agent of electroweak symmetry breaking remains elusive. The simplest model of symmetry breaking, called 
the Minimal Standard Model, uses a single complex doublet of fundamental scalars and is the benchmark for 
studies of the symmetry breaking sector of the theory. Direct search limits from LEP tell us that the Higgs 
mass is greater than ~ 114 GeV. Fits to precision electroweak data prefer a mass well below the direct search 
limit although the 95% confidence level upper limit is somewhat greater than 200 GeV. 

Higgs boson production at hadron colliders is dominated by the gluon fusion mechanism. However, experi- 
ments must not only produce Higgs bosons, they must also detect them. With a center-of-mass energy of 2 TeV, 
the Fermilab Tevatron is primarily sensitive to a Higgs boson with mass below the threshold for decay into W 
boson pairs. In this case, the Higgs will decay almost exclusively into bb pairs which will be undetectable on 
top of an enormous QCD background. Since the total cross section is too small to permit the use of rare decay 
modes, a light Higgs can only be detected through associated production with a W or Z boson. Only if the 
Higgs is sufficiently massive that the WW* channel begins to open up, will inclusive production via the gluon 
fusion mechanism be useful in the Tevatron Higgs search. 

At the CERN LHC, however, gluon fusion will be the discovery channel for the Higgs. The cross section for 
light Higgs boson production will be sufficiently large that the rare decay H — > 77 can be used up to the point 
that the WW* channel begins to open up. From that point on, the diboson decays provide a very robust signal. 

II. METHODS 

The Higgs boson couples to mass, which presents a problem for hadronic production. Gluons are massless 
and therefore do not couple directly to the Higgs at all, while the quarks that make up the proton have very 
tiny masses. Therefore, the dominant production mechanism is gluon fusion via virtual top quark loops. In 
the limit that the top quark is very heavy, we can integrate out the top and formulate an effective Lagrangian 
coupling the Higgs boson to the light quarks and gluons jl], |2[ If we take the light quarks to be massless, the 
effective Lagrangian takes the form 

Cos = Ctia^HG^G^, (1) 

where G" is the gluon field strength tensor. The coefficient function C\(a s ) has been computed to order Q. 

The use of the effective Lagrangian allows us to replace massive loop diagrams with point-like interactions. 
Next-to-leading order (NLO) corrections to inclusive Higgs production have been computed using both the 
effective Lagrangian || || and the full theory (t], ||. One expects that the effective Lagrangian will work very 
well if the Higgs mass is much smaller than twice the top mass but that it will be unreliable for larger masses. 
In fact, it was found that at NLO the effective Lagrangian does indeed agree very well with the full calculation 
below the top threshold and was even found to agree to within 10% for Higgs masses as large as 1 TeV. 
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FIG. 1: Representative diagrams of virtual corrections, single real radiation and double real radiation. 



It was also found that the NLO corrections are very large, of order 70 — 100%. Such large corrections clearly 
call for the evaluation of still higher-order terms in order to arrive at a solid theoretical understanding of the 
process. Since the effective Lagrangian seems to be a valid approximation, especially in the phenomenologically 
interesting region of Higgs boson masses below 200 GeV, we have embarked on an effort to compute the full 
next-to-next-to-leading order (NNLO) corrections in the heavy top limit. In this talk, we will present results for 
soft plus virtual corrections to inclusive Higgs production M, |l^, ^l) . These terms are not expected to dominate 
the full result and for this reason we also discuss an approximation of the formally sub-leading but numerically 
dominant contribution 1112] . 



III. THE SOFT APPROXIMATION AND BEYOND 



There are three distinct contributions to inclusive Higgs production at NNLO (see Figure |l|): Virtual correc- 
tions to two loops, single real radiation to one loop and double real radiation at tree level. These three channels 
produce radiative corrections that fall into three categories, depending on their functional dependence on the 
fraction x = M^/s of the center-of-mass energy squared of the scattering process that goes into producing the 
Higgs boson. 
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In the virtual corrections, all of the energy goes into Higgs boson production, so these terms contribute only to 
the 8(1 — x) correction. Real emission processes generate terms like (1 — x) n_m£ where n > — 1 is an integer 
and e is the dimensional regularization parameter where space-time is taken to be D = 4 — 2e dimensions. 
These processes contribute to the a and b n coefficients in equation (|^) by expanding terms like (1 — x) _1_me as 
distributions 
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In the soft limit, there would be no energy carried away by real emission and only the <5(1 — x) term would 



be kept. However, the 



ln"(l-z) 



terms are directly connected to the <5(1 — x) terms through equation (0) and 



in canceling the infrared poles proportional to 6(1 — x) we get these terms for free so they are kept as part of 
the soft approximation. 

While the soft approximation keeps the formally leading terms, it was found that at NLO it is a poor 
approximation. It is actually the sub-leading c n (n = 0,1 at NLO) terms that dominate the cross section. 
At NNLO, the c„ terms are again expected to dominate. Kramer, Laenen and Spira have used collinear 
resummation to derive approximate NNLO results for a,b n and c„. We expect their resummation to give the 
correct values for the coefficients 63, 62 and C3. The other coefficients require additional calculation, higher order 
resummation coefficients or, for the remaining c„, receive non-collinear contributions and we do not expect the 
approximation to be accurate. For the a and b n terms which we have computed directly, these expectations are 
fulfilled, giving us confidence that the dominant term, C3 is indeed accurate. 

This gives us a range of possibilities for estimating the full NNLO correction. In Figure || we show three 
approximations in addition to the soft limit: 

1) Use C3, C2, ci and cq from Ref. fMf , 

2) Use C3 from Ref. [12 and generate sub-leading ln"(l — x) terms by expanding xb n — » b n + (1 — x)b n , 

3) Use C3 from Ref. [12 and drop all sub-leading ln™(l — x) terms. 
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NNLO K-factor (14TeV) [CTEQ5] 



NNLO K-factor (2 TeV) [CTEQ5] 
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FIG. 2: NNLO K-factor as a function of Higgs boson mass at (a) s/l = 14 TeV and (b) s/I = 2 TeV. The solid line 
corresponds to approximation 1 above, the dashed line to approximation 2 and the dash-dot line to approximation 3. 
The dotted line represents the soft approximation. 



Note that in order to truly estimate the NNLO cross section, one needs NNLO parton distribution functions 
(PDFs). Unfortunately, the necessary ingredients for producing NNLO PDFs are still being developed. Ap- 
proximate NNLO PDFs have been produced, but at the time of this work they are not yet publicly available. 
We therefore use the CTEQ5 NLO parton distribution functions [[l3) and acknowledge the inconsistency. 

There are two outstanding features of Figure |[ the formally sub-leading ln 3 (l — x) terms dominate the 
corrections, and the size of the corrections is very large. One expects that using NNLO PDFs will reduce the 
magnitude of the correction by ~ 10%, but it will still be very large. We can take the spread between these 
approximations as an estimate of the uncertainty due to the uncalculated terms. 



IV. CONCLUSIONS 



We have described a calculation of the soft plus virtual NNLO corrections to inclusive Higgs production and 
estimates of the full NNLO correction based on collinear resummation. While the soft plus virtual terms are 
perturbatively well-behaved, the leading non-soft terms dominate the cross section and give rise to very large 
corrections. At this time, the two most important questions concerning inclusive Higgs production are 1) What 
is the precise value of the NNLO K-factor? and 2) How reliable is the NNLO K-factor with respect to even 
higher order corrections? The first question can be answered by completing the full NNLO calculation and 



this work is underway 
properties of the Higgs 



14j. The second question, which is crucial for determining the precision to which the 
Doson can be measured at the LHC, requires further investigation. 



[1] A. I. Vainshtein, M. B. Voloshin, V. I. Zakharov, and M. A. Shifman, Yad. Fiz. 30, 1368 (1979), [Sov. J. Nucl. Phys. 
30, 711 (1979)]. 

[2] A. I. Vainshtein, V. I. Zakharov, and M. A. Shifman, Usp. Fiz. Nauk 131, 537 (1980), [Sov. Phys. Usp. 23, 429 
(1980)]. 

[3] M. B. Voloshin, Yad. Fiz. 44, 738 (1986), [Sov. J. Nucl. Phys. 44, 478 (1986)]. 

[4] K. G. Chetyrkin, B. A. Kniehl, and M. Steinhauser, Phys. Rev. Lett. 79, 353 (1997), hep-ph/9705240. 
[5] S. Dawson, Nucl. Phys. B359, 283 (1991). 

[6] A. Djouadi, M. Spira, and P. M. Zerwas, Phys. Lett. B264, 440 (1991). 

[7] D. Graudenz, M. Spira, and P. M. Zerwas, Phys. Rev. Lett. 70, 1372 (1993). 

[8] M. Spira, A. Djouadi, D. Graudenz, and P. M. Zerwas, Nucl. Phys. B453, 17 (1995), hep-ph/9504378. 

[9] R. V. Harlander, Phys. Lett. B492, 74 (2000), hep-ph/0007289. 
[10] R. V. Harlander and W. B. Kilgore, Phys. Rev. D64, 013015 (2001), hep-ph/0102241. 
[11] S. Catani, D. de Florian, and M. Grazzini, JHEP 05, 025 (2001), hep-ph/0102227. 
[12] M. Kramer, E. Laenen, and M. Spira, Nucl. Phys. B511, 523 (1998), hep-ph/9611272. 
[13] H. L. Lai et al. (CTEQ), Eur. Phys. J. C12, 375 (2000), hep-ph/9903282. 
[14] R. V. Harlander and W. B. Kilgore (2001), in Preparation. 



